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Abstract 

The weak form factor for B — » K\b where K\b is the l P\ state is calculated in the light-cone 
sum rules (LCSR). Combining the quark model result for the form factor of B — » K\a with K\a 
being the 3 Pi state, we have larger values for B — > K\ form factors than the previous LCSR results. 
The increased form factors reduce the discrepancy between theory and the experimental data for 
B — > if 17. Some phenomenological meanings are also discussed. 

PACS numbers: 12.38.Lg, 13.20.He 



*Electronic address: jplee@kias.re.kr 

1 



I. INTRODUCTION 



Radiative B decays to K are a rich laboratory for the standard model and new physics. 
Especially, B — > K*j is well understood theoretically via b — > sj transition as well as 
experimentally. Recently, higher resonant kaons are observed by CLEO and B factories 
For example, BELLE collaboration has measured the radiative B — > K\ decays for the first 
time 1^: 

B(B + -> if+(1270)7) = (4-28 ± 0.94 ± 0.43) x 10" 5 , (1) 
B(B + -> #+(1400)7) < 1-44 x 10~ 5 , 

where Ki is the orbitally excited axial vector meson. In the theoretical side, recent de- 
velopments of the QCD factorization (QCDF) j^] makes it possible to calculate the hard 
spectator contributions systematically in a factorized form through the convolution at the 



heavy quark limit. B —> K ' 7 is already studied in this line jj, [5, 0, One good point 
about K\ is that there are lots of things shared with B —>■ K* / y. Basically both of them are 
governed by b — > sj. And the distribution amplitudes (DA) of K* and Ki are same except 
the overall factor of 75 which makes few differences in many calculations. 

A straightforward extension of the analysis for B — > K*^ to B — > Kij was given in 
. But the BELLE measurements of Eq. © reveal that theory predicts much smaller 
branching ratio than data |9J, [10J ■ This is an opposite situation to that of B — > K*^ where 
theory predicts larger branching ratio. Considering the resemblance between K* and K\, 
it is quite unlikely that the same theoretical framework would produce discrepancies with 
experiment in a reversed way. 

In the previous analysis the main uncertainty of theory lies in the nonperturbative form 



factors. Ref. |8 



factors 



In 



relies on the light-cone sum rule (LCSR) results for the B — > Ki form 
111 ] only the leading twist DAs are considered without any non-asymptotic 



contributions. In this paper we revisit the B — > K\ form factors in the LCSR. There are 



three improvements compared to ll|. First, higher twist DAs are included; second, non- 
asymptotic contributions are also considered; third, terms proportional to m\, where tjia is 
the mass of axial meson, are not neglected. 



For B — > K* form factors, the LCSR results are updated [12j up to the one-loop correc 



tions to twist-2,3 contributions and leading order twist-4. It is thus legitimate to improve 
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the theoretical accuracy for B — > K\ form factors. 

It is believed that the physical If i (1270) and i^i(1400) states are the mixtures of angular 
momentum eigenstates l P\ (Kib) and 3 Pi (Kia)- The mixing angle is not known precisely, 
but is close to the maximal. This is a very natural and convenient way to explain the 
suppression of one decay mode compared to the other. For the suggestive angles 6 = 
±37° . ±5 8° [l3| . negative ones are disfavored by (J2J. 

In some of the Gegenbauer moments of K\ B DAs are calculated by LCSR. With this 
information, we explicitly calculate the B — > K\b form factor in LCSR. Since K\b and K\a 
have different G-parity, their Gegenbauer expansion will not be the same. Future study on 
Kia is necessary to reinforce the reliability of current work. We use the results from model 
calculations for K\a to give B — > K\ form factors. This form factor will also be available 
for nonleptonic decay modes Q] 

The paper is organized as follows. In the next section the weak form factors and axial 
vector meson DAs are defined. The LCSR evaluation is given in Sec. III. Section IV deals 
with the LCSR results. In Sec. V, some discussions about the results and their meanings 
appear. Conclusions are also added at the end of this section. 



II. FORM FACTORS AND DISTRIBUTION AMPLITUDES 

For the axial vector A(pA,e), where pa (e) is the momentum (polarization) of A, the 



relevant B — > A transition matrix elements are defined as 



(A(p A ,e)\qia^b\B{p B )) (2) 
= F A (q") [(e* • q)(p B +p A ), ~ <(m| - m\)] + F_V) [(e* • q)q, - etf] 

+^^-(e*-q)[(m 2 B-m A )q,-(pB + PA)q 2 ] , 

VTLBlTiA 

(A(pA,e)\qia, ul5 q»b\B(pB)) = -iF A {q 2 )e^e w q a {p A + p B f , (3) 



where q = Pb — Pa and 777.5(7??.^) is the B (axial vector) meson mass. We use e i23 = +1- 
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The distribution amplitudes (DA) of the axial vector meson are given by 



14, 



17 



(A(p A ,e)\qi(y)-f^ 5 q2{x)\0) 
= ifAmA [ du e i{ - upy+Upx) 



(4) 



e • z 



p ■ z~ 



-Pu 



e • z 



- ( f> ll (u)+2g i l\u)-g 3 (u) 



(A(p A , e)\q x (y)(T^ 5 q 2 (x)\0) 
Jo 



(6) 



e • z 



(p ■ z)' 



™\{p A uZv-pt z p) 



lcf> ± (u) + hf(u)-±h 3 (u) 



m 



A ( * 



2p ■ z 



z v - e*,z M ) [-(f>x(u) + h 3 (u)} 



(A(p A , e)\qx(y) l5 q,(x)\0) = f^m 2 A (e* • z) j du ^\u) 
Here z = y — x and 



P^ =Pa 



(7) 



(8) 



2p A ■ z 

is the light-like vector, and u — 1 — u. The DAs 0y (twist-2), gy , gy (twist-3), and g 3 
(twist-4) are anti-symmetric under the change u — > u while (j)± (twist-2), h$, h { ° ] (twist-3), 
and h 3 (twist-4) are symmetric in SU(3) limit, because of the G-parity. Thus 



Jo 



duf{u) = , for / 



(a) {v) 

9± , 9± , 93- 



(9) 



The leading twist DAs are expanded with the Gegenbauer polynomials. In general, we can 
expand 



I u 



6uu aj C\ 



3/2 



U — U) 



(10) 



1=0 



(f>±(u) = 6uu 



l + ^Cf /2 (u-u) 



1=0 
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For twist-3 DAs, the Wandzura-Wilczek type approximation will be used; 

dv 



ffVW 

( a ) I \ 
9±(u) 



u I dv 









dv 


V Ju 










— h U 


r 




a 



h 



(t) 



(u — u) 



is) 



(ii) 



h 



The twist-4 DAs will not be considered afterwards. The first few Gegenbauer coefficients 
are recently calculated by QCD sum rules (l4|. 



III. SUM RULE EVALUATION 

The main point of LCSR is to evaluate the two point correlation function: 

II4 = if d A x e- ip *- x (A(p A ,e)\T [j(0)jt (s)] |0) . (12) 

Here j B (x) = 6(x)i7 5 g(x) is the interpolating current for B meson, and J(y) = q{y)Tb{y) is 
the heavy-to-light current with V being an appropriate gamma matrices. To establish the 
sum rule, one calculates 11^ in two ways. On one hand, 11^ is described in terms of hadronic 
observables. We call this n^ ad . Explicitly, 

(A\J(0)\B)(B\jU0)\0) 



n 



had 



m% ~p 2 B + ie 



res. 



(13) 



where the first term is the B meson contribution and (res.) is the higher resonance one. 
The term (A\J(0)\B) defines the transition form factor while 



m. 



(B\U l5 q\Q) = -£f B , 



(14) 



is proportional to the B meson decay constant fs- Here n^ ad is considered as an analytic 
function of p B . Using the dispersion relation, 



n 



had 



00 p had (s) 

ds — 



1 ' 



(15) 
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where p had (s) is the spectral density function. This is another expression of Eq. (|13jh from 
which we can extract the form of p had . 

On the other hand, 11^ can be written by quarks and gluons, and hence by light-cone 
distribution amplitudes (LCDAs). We call this I1 LC . From the dispersion relation, 

Jm\ 5 Pb 

i r° ImlT^s) 

CLS 



2 

71 J ml s ~ Pb 



where the imaginary part of Il^ c will be expressed by the LCDAs. At this stage, one assumes 
the quark-hadron duality for (res.) in Eq. (JT3J) as 

res. = -ds / A -^ , 17 

J S0 s- p B 

up to possible subtractions. Here so is the continuum threshold from which higher multi- 
particle states begin. In the numerical analysis, s is considered as a hadronic parameter. 
Combining all this, one arrives at 

(A\J(0)\B)(B\j B (0)\0) _ 1 [ S0 ds Imny (g ) ^ 
m\-p\ 7T J m 2 s-p 2 B 

After the Borel transformation over p B , we have the final expression for the sum rule: 

i rso 

e- m ^ T (A\J(0)\B)(B\j B (0)\0) = - ds e^ T lmU L A c (s) , (19) 

^ J ml 

where T is the Borel parameter. 

Among the three form factors F A Q (q 2 ), the most important one is F A (q 2 = 0) since only it 
is responsible for the radiative decay of B — ► K\. Also, it can be shown that F A (0) = F A (0) 
[l^ . To extract F A , we find it convenient to choose J(0) = qicr^^b. The left- hand- side 
(L.H.S.) of Eq. (HU) is simply 



(L.H.S.) = -iF A {q 2 )^ a ^q a {p B +p A f (^fn) e~"^ T 



(20) 



The right-hand-side (R.H.S.) of Eq. ()19j) is rather involved. After contracting the bb quarks, 

(R.H.S.) = - ds e~ s l T Im / d 4 x / -f = (21) 

vr J ml J J (2vr) 4 k 2 - m 2 + ze 

{A\q{<S)a^q v 1tq{x)\$) + rn h {A\q{<S)a^q v q{x)\Qi) 
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The two matrix elements in the above equation can be written, after some gamma matrix 
algebra, in terms of the LCDAs, Eqs. (jU)-©. In Eqs. (jU)-©, the position coordinate x can 
be replaced effectively by 







id{up)^ ' 

which is guaranteed by the presence of e tupx . On the other hand, for the factor 1/p ■ x 
1 , 



p ■ x 
1 



( // ; - - i j dv (j){v) , 
'o 

(u) — > i 2 / dv dw 4>(w) 
Jo Jo 



(22) 



(23) 



(a) (v) \ 

g± g±, 93) 



(p ■ x) 2 

where the surface terms are vanishing. In this way, one can remove x-dependence in (R.H.S.) 
except in the exponent. Thus the integration over x yields a delta function, ~ 5 4 (k— ps+up). 
Another delta function appears in the imaginary part of l/{k 2 — m 2 B + ie). Combining all 
together, one arrives at 



(R.H.S.) 



(24) 



«W f *VPi / ds e s/T [ du\ f A m A \g { l\u) - 25 s + (-us + um\ - (1 + u)q 2 )5' s 



f A m A $||(m)5 s + M5fJ ) (M)5 s -m 2 A G 3 (u)8' s + f A m b <f>±(u)6 s - m 2 A H 3 (u)S' s 



Here we use the short-hand notation, 8 S = 5(s — m 2 — 2up ■ Pb), and the differentiation is 
with respect to up. It is understood that at the final stage of calculation, p — > p A . And the 
newly defined functions are 



$|| (u) 
G 3 {u) 
H 3 (u) 



dv 



o 



M v ) ~9 { l\ v 



(25) 



dv / dw 
ii Jo 



dv 



4> ll (w) + 2g^(w)-g 3 {w) 
h 3 (v) - (f)±(v) 



Equating Eqs. (|20|) and (|25|h after a little algebra, we have the final expression for the form 
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factor F+(q 2 

fHo) 



where 



hadronic information (in GeV) 


Gegenbauer moments (at 1 GeV) 


m B 5.27 


all 


0.26 


m b 4.8 


1 


-1.75 


fs 0.161 


a| 


0.13 


m A 1.370 




-0.13 


/a 0.195 




-0.02 




a 3 


-0.02 



TABLE I: Input values. Gegenbauer moments are from 



0) 

1 



e m %/ T 



m b 



m 2 B l 



B 



fAtriA 



so/T 



s + m A 



s - m A j a) 
4 



9 ^(uo)+mf 3{Uo) 



u 



+fAm A I — exp 



"0 



U 



m\ + um\ 
vT 



x 



uT + m\ + (1 — 2u)m 



A (a) 



-f A m b e 



AuT 

2 

-sp/T m A 



g^(u)-^(u)~ug ( l'(u) + 



(«)/ \ . m A G 3 (u) 



a 



H 3 (u ) 



s + m A u 
+f A m b / — exp 



u 



m 2 + um A 
vT 



<f> ± (u) 



m 2 A H 3 (u) 
T u 



Uq = 



m b + m \ 

So + m\ 



(26) 



(27) 



IV. RESULTS 



In what follows, only the case where q 2 = is considered. The basic input constants 
are summarized in Table HV1 The LCSR involves two important parameters, the continuum 
threshold So and the Borel parameter T. Naively thinking, the continuum threshold is 
roughly 

s ~ (m B + A) 2 = (2m B - m b ) 2 , (28) 
where A = m B — m b . Numerically, 



s ~ (2m B - m b ) 2 = s* » 33 GeV 2 
8 



(29) 



for rriB = 5.27 GeV and m& = 4.8 GeV is consistent with literatures Q|. We take this value 
as a starting point to fix sq. 

In principle, is independent of the unphysical Borel parameter T. But in reality there 
is a sum rule window of T where a physical quantity is stable. If T is too small, then the 
higher twist terms proportional to 1/T n (n = 1,2, • • •) become too large. One requires, for 
example, 

(^r terms) _ . s 

This condition imposes the lower bound of T. The number 30% might be changed, but we 
adopt this value here. On the other hand, if T is too large, then the contributions from the 
continuum states become too large. We require that 

i ds e- s / T Imn(s) 

tfr^™ 5 30% ' (31) 

This constraint imposes the upper bound of T. Note that the condition of Eq. (|3*T|) is used 



m 



12j to determine the lower bound of continuum threshold, Sq. In this analysis, however, 
we start with sq = to determine the sum rule window, and then we fix so from the best 
stability of within the sum rule window. 

From Eqs. (pffij) and (J3T|) with sq = s*, we have 

6.8 < T < 21.7 (GeV 2 ) . (32) 

This window has overlaps with that of ^l|, but not with that of [3] where only vector 
mesons are considered. As an illustration, plots of over T for various Sq around s* are 
given in Fig. ^ To find the best value of sq, we impose a simple condition. We scan so which 
minimize the value F^(T C + 5 GeV 2 ) - F^(T C - 5 GeV 2 ), where T c is the central value of T 
within the sum rule window. We find that the best value of Sq is 

s b = 34.3 GeV 2 . (33) 

Plots of Fjj 4 for various sq around Sb are shown in Fig. 121 A closer look of Fig. 121 is given in 
Fig. and 3-dimensional plot of F+ against Sq and T is given in Fig. 0] From the above 
analysis, we get 

i^(g 2 = 0;T c ) = 0.256™°, (34) 
where the errors are from the variation of Sq around s& by ±1 GeV 2 . 
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FIG. 1: F$(0) vs T for various sq around 




s b + 2 (GeV 2 ) 

_s 6 - 1 
s 6 -2 



7.5 



10 12.5 15 17.5 20 22.5 



T (GeV 2 ) 



FIG. 2: F£ (0) vs T for various so around s&. 

The observed axial kaons i£i(1270) and (1400) are mixtures of X P\ and 3 Pi states. 
Their form factors are related via mixing angle 9 as |l3j, [l9( 

^^(1270) = ^3 sin + F A CQsei (35) 

where i — 0, +, — . Here, F/ 43 are the 3 Pi form factors. We use the result of F+ 3 (q 2 = 
0) = 0.11. The mixing angle 9 is not yet fixed precisely. Ref. [l^] suggests 9 = ±37°, ±58°. 
Table HP shows the values of p^ Kl ( 127QS > anc [ pB-^Ki(uoo) ^ these angles. For negative 
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FIG. 3: A closer look of Fig. El 




FIG. 4: 3-dimensional plot of F+. 
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—B-vKifmO) 



F 



73^^1(1400) 
+ 



37° 



-37° 58° 



-58° 



0.271 0.138 0.229 0.042 
-0.066 0.242 -0.159 0.276 



TABLE II: p^^Om and F ^ l( i400) for 
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+(3) 


+(3)+(2) 


+(3)+(2)+(l) 


0.187 


0.237 


0.155 


0.256 



TABLE III: Contributions of new improvements. The parameter set used here is the same as that 
in the previous section. Conditions (1), (2), and (3) are explained in the text. 

angles, we find 

F B->Ki(1270) < ^B^Xi(1400) ^ 

Since other parameters of the branching ratio are not so different in B — > (1270)7 an d 
B -> ^(1400)7, one expects B[B -> 7^(1270)7) < B{B -> (1400)7) for the negative 
mixing angles. This is not consistent with the experimental data. 

V. DISCUSSIONS AND CONCLUSIONS 

As discussed in the discrepancy between theory and experiment for B{B — > 
i^i(1270, 1400)7 is mainly due to the smallness of the relevant form factors. If there is 
no mixing (i.e., 6 = 0), then pf^ Kl ( 1270 ) = pA _ q 256. This is considerably larger than 
the previous LCSR result of (llj, F^J^^ 1270 ^ = 0.14 ± 0.03. The mixing effects are only 
5.7% and —10.6% for 9 = 37°, 58°, respectively. In jll|, only the asymptotic form of leading 
twist DA, 

07 (u) = 6u(l - u) , (37) 



ml + um\ 



contributes to the sum rule. According to Eq. (31) of [11 



nT 0l Sy (^) • (38) 

It should be compared with Eq. (|27j). Eq. ()27|1 improves Eq. ()38|) in three ways: (1) higher 
twist DAs are included; (2) non-asymptotic contributions are also included; (3) there is no 
term proportional to m\ in g afir (0). With the parameter set used in the previous section, 
we have Safir (0; s = s&; T = T c ) = 0.187. This is lager than the value of F^J^^ 1270 ^ = 
0.14 ± 0.03. But if we take the sum rule window of Borel parameter adopted in jll|. 
F± g a g r (0; s = Sb;T = 7.5 GeV 2 ) = 0.151, which assures the consistency of the present 
analysis. We can check how much the new improvements contribute to the form factor. The 
results are summarized in Table IIII1 One finds that non-asymptotic and higher-twist DA 
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e 


0° 


37° 


58° 


experiment [2] 


B(B° - 


-> iTi(1270)°7) 


2.81 


3.14 


2.27 


4.28 


B(B° - 


-» Ki(1400)°7) 


0.52 


0.14 


0.91 


< 1.44 



TABLE IV: Branching ratios for various mixing angles in units of 10 5 . 

contributions as well as non-zero mass effects are considerable. 

The increase of the form factor will reduce the discrepancy between the theoretical pre- 
dictions and experimental data Q|. At next-to-leading order of a s , the branching ratio of 



B — > Ki'j is given by |8, 



1 F^(0)(-0.385-z0.014) + (/|/GeV)(-0.024-i0.022) 

m B J 

(39) 

The resulting branching ratios are given in Table IIVI The enhancement is significant and 
the theoretical prediction becomes closer to the experimental data compared to the previous 
analysis [8], though there is still a gap. 

There are a few possibilities to improve further. Firstly, the precisions are different 
between Eqs. (|39|) and (|2Tj) . Eq. (|39|) contains the hard spectator interactions which appear 
as a convolution between the jet function and the meson DAs. The DAs contributing to 
Eq. ()39|) are leading twist ones and of asymptotic form. It is thus necessary to include 
higher twist and non-asymptotic contributions in Eq. (}39|) at the same accuracy as was 
done in this work. Also, Eq. (|2"Tj) contains terms proportional to m\, but Eq. (|3!?|) is the 
result of heavy quark limit. One can easily expect that the next-to-leading order (NLO) 
of Aqcd/^& corrections to the QCDF framework might include the terms of m^/mj,, but 
there is no systematics so far. The hard spectator interactions are given by the convolution 
of hard kernel and meson DAs. Similar non-zero m 2 A terms will also appear in the axial 
vector DAs to affect the hard spectator interactions. But this effect is not expected to be 
large. According to j8j, the hard spectator contributions amount to roughly about 5% at 
the amplitude level. 

Secondly, F+ 3 could be larger. Actually there is no clue about the size of but it 
might be that Fjj 43 is comparable in size to F+, just as in ^J. If this is the case, then the 
form factor can be enhanced via mixing 

f b-*Kx{wo) ^ Q 256 x ( sin6 + cos ^ = o_35 _ 0M ) ( 40 ) 

13 



for 6 = 37°, 58°, which results in a large branching ratio, 

B{B° -> if 1 (1270) 7) « 5.3 x 1(T 5 . (41) 

In conclusion, we have calculated B — > form factor in LCSR. This analysis improves 
the previous one in a few respects by including higher twist DAs, non-asymptotic contribu- 
tions, and non-zero m 2 A terms. The value is rather larger than the previous calculation and 
that from the quark model result. Larger value is well accommodated to the experimental 
data. One needs more information about K\a and the mixing angle to reduce theoretical 
uncertainties. To go beyond the current work, one can include the NLO of a s which might 
not be so different from that for B — > K* [12]. And the study of B — > Ki"/ at higher 
accuracy comparable to this work will be necessary. 
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